arXiv: 1502.00001v3 [astro-ph.SR] 28Jul2015 


000, 000-000 (0000) Printed 30 July 2015 (MN style file v2.2) 


loss’s First Supernova? New Limits on the “Impostor” 
SN 1997bs 


Scott M. Adams^ & C.S. Kochanek^’^ 

^ Deptartment of Astronomy, The Ohio State University, IfO W. 18th Ave., Columbus, OH 43210, USA 

^ Center for Cosmology and AstroParticle Physics (CCAPP), The Ohio State University, 191 W. Woodruff Ave., Columbus, OH 43210, USA 
E-mail: sadams@astronomy.ohio-state.edu 


30 July 2015 


ABSTRACT 

We present new, late-time Hubble Space Telescope and Spitzer Space Telescope ob¬ 
servations of the archetypal supernova (SN) impostor SN 1997bs. We show that SN 
1997bs remains much fainter than its progenitor, posing a challenge for the canonical 
picture of late-time obscuration by dust forming in a shell ejected during the tran¬ 
sient. The possibility that the star survived cloaked behind a dusty, steady wind is 
also disfavored. The simplest explanation is that SN 1997bs was a subluminous Type 
Iln SN, although it is currently impossible to rule out the possibility that the star sur¬ 
vived either behind an obscuring dusty shell > IMq or with a significantly decreased 
intrinsic luminosity. Continued monitoring of the system is the most promising means 
of ultimately resolving this question. 

Key words: stars: evolution - supergiants - supernovae: general - supernovae: indi¬ 
vidual: SN 1997bs 


1 INTRODUCTION 

Supernova (SN) impostors (Van Dyk et al. 2000) are a 
class of stellar transients characterized by Type Iln spectra 
(narrow hydrogen emission lines) with lower peak luminosi¬ 
ties {Mv — —13) than typical core collapse SNe (ccSNe; 
Mv — —17). The class is heterogeneous, including luminous 
variable stars (e.g. SN 2002kg; Weis & Romans 2005) and 
the SN 2008S class of transients (Prieto et al. 2008; Thomp¬ 
son et al. 2009; Kochanek 2011a), but a subset of the events 
are generally connected to some type of eruptive transient 
associated with the phenomenology of luminous blue vari¬ 
able stars (LBVs; Humphreys & Davidson 1994). Still, the 
true nature of SN impostors is debated, including whether 
they are non-terminal eruptions or actual SNe (see Smith 
et al. 2011; Kochanek et al. 2012). 

The rate of SN impostors attributed to LBV eruptions is 
likely ~ 20%-60% of the ccSN rate (Thompson et al. 2009). 
Thus, the true nature of these transients could have impor¬ 
tant consequences for the rate of SNe. For example, the SN 
rate appears to be less than the massive-star formation rate 
(Horiuchi et al. 2011). There are several possible solutions to 
this mismatch. Non-local SN surveys could be significantly 
incomplete (e.g., Botticella et al. 2012) or there could be a 
significant rate of failed SNe (Kochanek et al. 2008; Gerke 
et al. 2015). The latter solution could also help to explain 
the lack of higher mass SN progenitors (> 15 M©; Kochanek 


et al. 2008; Smartt 2009) and the black hole mass function 
(Kochanek 2014c, 2015). The final option is that some of 
these lower luminosity transients are in fact SNe (Horiuchi 
et al. 2011). 

On the other hand, if SN impostors are non-terminal 
events, they may be the dominant mode of mass-loss in mas¬ 
sive stars (Smith & Owocki 2006). This would be important 
because eruptive mass-loss is unaccounted for in stellar evo¬ 
lution models and little is known about the mechanism, duty 
cycle, mass dependence, or total mass-loss in such events. 
Recently Khan et al. (2013, 2015) surveyed nearby galax¬ 
ies and found that luminous stars encased in ejected dusty 
shells analogous to rj Carinae are not common enough for 
eruptions to be the dominant mass-loss channel, represent¬ 
ing ~ 10% of mass-loss rather than ~ 50%. 

Inspired by the phenomenology of rj Carinae, the most- 
widely accepted picture of SN impostors is that these events 
are the result of non-terminal eruptions that eject a signif¬ 
icant amount of mass that then may obscure the surviving 
star (see e.g., Humphreys et al. 1999). The phenomenology 
expected for a dusty shell is fairly simple - the optical depth 
peaks as the shell passes through the dust formation radius 
and then declines as ~ with time (Kochanek 2011b). 
Kochanek et al. (2012) examined several SN impostors with 
archival Hubble Space Telescope (HST), Spitzer Space Tele- 
seope (SST), and Large Binocular Telescope (LBT) data 
and found that none of the transients (for which there are 


0000 



2 Adams & Kochanek 


sufficient data to reach a conclusion) are consistent with the 
shell ejection scenario. If they are simply outbursts, the data 
are more consistent with a scenario where the transient is a 
sign post that the star is transitioning from a low to a high 
mass-loss rate, resulting in a longer lived dusty wind. Still, 
nothing in the available data ruled out the possibility that 
some of these systems are not impostors at all, but instead 
are true SNe. 

The basic problem is that SN impostors have been lit¬ 
tle studied after the brief optical transient even though it 
is really the late-time phenomenology that is the key to un¬ 
derstanding these events. For example, it would be useful 
to simply establish the continued presence of a star. In this 
paper we utilize new late-time observations to investigate 
the nature of the archetypal impostor SN 1997bs. 

SN 1997bs, discovered on 1997 April 15, was the first SN 
reported by the Lick Observatory Supernova Search (LOSS; 
Treffers et al. 1997). It was found in NGC 3627 and was 
classified as a Type Iln SN (the spectrum was dominated 
by Balmer emission lines with full width at half-maximum 
~ 1000 kms“^). Van Dyk et al. (1999) identified a candi¬ 
date progenitor with mpeoew — 22.86 ± 0.16 mag, which 
corresponds to a luminosity of L* = lO'^'^-lO®'^ Lq for the 
temperature range T* = 7500-20,000 K. The event peaked 
at V ~ 17 and had faded to (7^55514^) V ~ 23.4 mag by 1998 
January 10. Over this period the transient became redder, 
evolving from V — I = 0.7 to V — I = 3 mag (Van Dyk 
et al. 2000). Van Dyk et al. (2000) suggested that the flat¬ 
tening of the late-time light curve at ~ 0.5 mag fainter than 
the progenitor was an indication that the star survived the 
explosion. 

In fact, it continued to fade, since Li et al. (2002) found 
that the SN was marginally detected in an HST /Wide Field 
Planetary Camera 2 (WFPC2) image taken in 2001, with 
F’hhhlT = 25.8 ± 0.3 mag (March 4) and that it was not 
detected at TSldlV to a limiting magnitude of about 25.0 
(February 24 and May 28). Li et al. (2002) posited that while 
the formation of dust in the ejecta could explain the contin¬ 
uous decline in the optical flux, this would be inconsistent 
with the color evolution. The SN would become progres¬ 
sively redder if dust was forming in the ejecta, but instead 
the source appeared to become bluer between early 1998 and 
early 2001. We present the full optical light curve in Fig. 1. 

Van Dyk & Matheson (2012) reported mid-infrared 
emission near the position of SN 1997bs in SST Infrared 
Array Camera (IRAC) images from 2004 May, with 32 /ijy 
at 3.6 /rm and 40 /rjy at 4.5 /rm, but no detections at 5.8 /rm 
and 8.0 /rm. Van Dyk & Matheson (2012) stated that this 
could be fitted by a T ~ 970 K blackbody with a radius 
R ~ 1.4 X 10^® cm and luminosity L ~ 3.1 x 10® Lq. 
This would correspond to an expansion speed of only ~ 
60 kms“^, far slower than the ~ 765 kms”'^ line width 
observed during the transient (Smith et al. 2011). Addition¬ 
ally, Van Dyk & Matheson (2012) found that the star is 
recovered at rriFsssw = 26.08 and m_F 8 i 4 W = 25.08 mag in 
i/Sr/Advanced Camera for Surveys (ACS) data obtained 
in 2009, although they do not report uncertainties. 

Kochanek et al. (2012) considered the same SST obser¬ 
vations but had a different interpretation. They found that 
most of the flux of the point source possibly present in the 
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Figure 1. Optical light curve of SN 1997bs. The black filled 
squares (red open triangles) show the evolution of the F-band (/- 
band) luminosity converted from the HST magnitudes reported 
by Van Dyk et al. (2000) and Li et al. (2002) as well as our mea¬ 
surements of archival HST images from 2004 and 2009. The last 
set of V and /-band luminosities represent the upper limits we 
find for the new HST data we present in this paper. The short 
lines extending from the LBT point (the open pentagon) represent 
the F-band variability and 3cr limits found by image subtraction 
of our LBT monitoring data. The horizontal dashed lines give the 
range of progenitor luminosities for T* = 7500 — 20, 000 K given 
the candidate progenitor detection of m^gogw — 22.86 ± 0.16 
mag reported by Van Dyk et al. (1999). 


4.5 /im image vanished in wavelength-differenced images^, 
suggesting that the source of the flux was not dusty. Fur¬ 
thermore, they also found that the evolution of the spectral 
energy distribution (SED) is inconsistent with the ejection 
of a single shell at the time of the transient. Keeping the star 
obscured at the time of the HST observations in 2001 would 
require a continued mass-loss rate of M ~ 10~^Mq yr~^. 
However, matching the candidate surviving star to the pro¬ 
genitor requires a temperature (T* > 10^ K), too high to 
allow the formation of dust (see Kochanek 2011a, 2014b). 
Kochanek et al. (2012) concludes that if the star found by 
Li et al. (2002) is the survivor of SN 1997bs, the most likely 
scenario is that the eruption ended before 2001 and the star 
would have to be unobscured in 2011. 

In this work, we revisit these questions with new HST, 
SST, and LBT observations. In §2 we present the data, ex¬ 
plain our identification of the source and its photometry, 
and present the methods we use to model and constrain 


^ Lfsing difference imaging methods to subtract the 3.6 and 
4.5 iim images removes all sources with the ‘Rayleigh-Jeans’ mid- 
IR spectral energy distributions of normal stars to leave only those 
with significant dust emission. Thus, the lack of a counterpart in 
the wavelength-differenced images indicates an absence of dust 
emission. 
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the existence of a surviving star. As can be seen from Fig. 
2 , no star is obviously visible, and there is certainly no 
source with a flux comparable to the progenitor. In §3 we 
present the results of modeling the SED of the source as 
obscured by an expanding shell or a steady-state wind and 
show that they are inconsistent with a surviving star as lu¬ 
minous as the progenitor. Finally in §4 we summarize our 
results and consider alternative explanations for the fate of 
SN 1997bs. When converting observables into physical quan¬ 
tities we adopt the Cepheid distance of 9.4 Mpc to NGC 
3627 from Freedman et al. (2001) and a Galactic extinction 
of E{B — V) = 0.04 from the Schlafiy & Finkbeiner (2011) 
recalibration of Schlegel et al. (1998). We note that this re¬ 
sults in a significantly smaller progenitor luminosity than 
that adopted by Van Dyk et al. (2000), who used an older 
Cepheid distance of 11.1 Mpc from Saha et al. (1999) and 
also included an estimate of galactic reddening from NGC 
3627 of E{B — V) = 0.21 mag. Though we do not include 
this galactic reddening, such reddening, to first order, will 
be taken into account by our circumstellar dust models (but 
see Kochanek et al. 2012). Additionally, most of our results 
are dependent primarily on the luminosity of the progeni¬ 
tor relative to a possible survivor, which is independent of 
foreground extinction and the adopted distance. 

2 DATA AND MODELS 
2.1 Data 

We utilize both new and archival HST data. For our pro¬ 
gram, we obtained new Wide Field Camera 3 (WFC3) UVIS 
F555IF and F814IF and IR FllQW and E160W images 
(GO-13477) taken in 2013 November. We also use public 
HST WFC3 UVIS F275W, F336W, FA38W, F555W, and 
FSldlU images taken for the Legacy ExtraGalactic UV Sur¬ 
vey (LEGUS; PI D. Calzetti, GO-13364) in 2014 Febru¬ 
ary and the multi-epoch archival WFPC2 ESSSIU images 
of NGC 3627 taken between late 1997 and early 1998 (PI 
A. Sandage, GO-6549) that contain the transient event to 
calibrate our astrometry and measure the position of SN 
1997bs. Additionally, we analyze archival ACS/Wide Field 
Channel (WFC) F555W and FSldlU images taken in 2009 
December (PI S. van Dyk, GO-11575) and 1^43511^ images 
taken in 2004 December (PI R. Chandar, GO-10402). 

We also make use of new and archival SST data. We co¬ 
added archival images of NGC 3627 (from program ID 159) 
taken in 2004 May and also co-added archival images taken 
in 2014 Feb-March-Aug (ID 10136) together with data from 
our program (ID 10001) taken in 2013 July using the mopex 
reduction package^. We also performed image subtraction of 
these data using ISIS (Alard & Lupton 1998; Alard 2000) to 
check for source variability. 

We have been monitoring NGC 3627 with the LET as 
part of a program searching for failed SNe (Kochanek et al. 
2008; Gerke et al. 2015). We use isis to check for optical 
source variability in the dozen epochs we have collected be¬ 
tween 2008 and 2014. 


^ http;//irsa.ipac.caltech.edu/data/SPITZER/docs/dataanal- 
ysistools / tools / mopex / 


2.2 Candidate Identification 

Proper alignment of all the data and a precise measure¬ 
ment of the coordinates of SN 1997bs are critical for correct 
candidate identification and photometry. We first aligned 
and stacked the archival drizzled HST WFPC2 F555IU im¬ 
ages using SExtractor (Berlin & Arnouts 1996), scamp 
(Berlin 2006), and SWARP (Berlin et al. 2002). We chose 
to register the astrometry of all our data to the drizzled 
HST WFG3 F’814IF image from 2013 November. We found 
the position of SN 1997bs in our reference WFC3 FSldlU 
frame by aligning the stacked archival WFPC2 F’555IU im¬ 
age with the GEOMAP task in iraf® using a matched coor¬ 
dinate list of several dozen sources within 20” of SN 1997bs. 
We estimate an uncertainty of 0/002 in the astrometry by 
measuring the rms in the SN position using different sub¬ 
samples of the matched coordinate list with different orders 
for the astrometric fits. We measure a centroid uncertainty 
of 0/004 based on the aperture photometry of SN 1997bs 
in the stacked archival frame. Adding these uncertainties in 
quadrature gives a total positional uncertainty of 0/004 (0.1 
pixels) in our reference image. 

We generate photometric catalogs for the HST WFG3 
data using the software package dolphot 2.0 (Dolphin 
2000)^. We largely used the same parameter settings as in 
Dalcanton et al. (2012), with the exceptions listed in Table 
1. We again used the drizzled HST WFC3 ESldlU image 
from 2013 November as the reference. Although dolphot 
has a built-in routine for aligning images, it was necessary 
to calculate initial pixel offsets using the tweakreg task in 
the DRIZZLEPAC package for the images taken on a different 
epoch (2014 February) than the reference. The rms align¬ 
ment residuals for the different images ranged from 0/003 
for the UVIS 1^55511^ and F81AW images to 0/03 for the 
UV images, and 0/16 for the IR images. We also performed 
aperture photometry with the PHOT task in iraf. The aper¬ 
ture corrections were calculated using the HST point spread 
function (PSF) models from tiny tim (Krist 1995; Krist 
et al. 2011)®. 

The closest source found by dolphot at /-band is 0/04 
from the expected position of SN 1997bs with a centroid un¬ 
certainty estimated from the reported S/N of the detection 
of 0/004. We also estimate the uncertainty by bootstrap re¬ 
sampling the set of six F81AW images that are co-added 
to form the reference image. The rms in the location of the 
closest source to the position of SN 1997bs using these dif¬ 
ferent reference images is 0/11. The rms of the source lo¬ 
cation found using dolphot with slightly different sets of 
parameters is ~ 0/03. Taken together, these uncertainties 
mean that the source position is consistent with the position 
of SN 1997bs. The photometry for the AGS/WFC images 
was found by again running dolphot with the same refer¬ 
ence image. We report the dolphot photometry of the HST 
data in Table 2, separately listing the statistical uncertain¬ 
ties from DOLPHOT and systematic uncertainties estimated 


^ IRAF is distributed by the National Optical Astronomy Observa¬ 
tory, which is operated by the Association of Universities for Re¬ 
search in Astronomy (AURA) under cooperative agreement with 
the National Science Foundation. 

http://americano.dolphinsim.com/dolphot/ 

® http://tinytim.stsci.edu/cgi-bin/tinytimweb.cgi 
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Table 1. DOLPHOT Parameters 


Description 

Parameter 

WFC3/UVIS 

WFC3/IR 

Photometry aperture radius 

RAper 

4 

3 

Inner sky radius 

RSkyO 

15 

8 

Outer sky radius 

RSkyl 

35 

20 

X-statistic aperture size 

RChi 

2 

1.5 

Sky Fit 

FitSky 

1 

1 

Spacing for sky measurement 

SkipSky 

2 

2 


Only parameters differing from those used by Dalcanton et al. (2012) are 
listed here. 



Figure 2. The region surrounding SN 1997bs in the new HST 
WFC3 images, the ACS F555VP image from 2009, and the pre-SN 
WFPC2 F555VP image. The radius of the green circles (O'.'l) is 
25 times our positional uncertainty and each image is ll'O across. 


from the standard deviations of the fluxes measured using a 
range of dolphot parameter settings. 

We also looked for an IR source at the position of SN 
1997bs. First we registered the SST images to the JIST 
data using the GEOMAP task in iraf with coordinates 
of matched point sonrces. The point sources were identi¬ 
fied by subtracting SST images convolved with a larger ker¬ 
nel from those convolved with a smaller kernel. The rms in 
the astrometric solution is ~ 0'.'12. There is no clear point 
source at the position of SN 1997bs in the SST data but 
we give upper limits based on aperture photometry in Ta¬ 
ble 2. Unfortunately SN 1997bs is coincident with diffuse 
IR emission from a spiral arm, making photometric mea¬ 
surements challenging. In order to minimize contamination 
from diffuse emission and other sources we use a 1'.'2 radius 
aperture with a l'.'2 — 3'.'6 radius sky annulus together with 
empirically-determined aperture corrections. While such a 
small sky annulus Includes significant point source flux, the 
background is non-uniform on larger scales. We also mea¬ 
sure the variability of SN 1997bs in the SST images taken 

000 


between 2004 and 2014 using image subtraction. While the 
3.6 fim flux stays constant, the 4.5 /rm flux declines by 24 
pjy. 

We use the LET data to place limits on optical variabil¬ 
ity at the coordinates of SN 1997bs. No source is detected 
at the position of SN 1997bs and with a dozen epochs span¬ 
ning 2008-2014 we find slopes consistent with zero in all 
filters (150 ± 420, -260 ± 450, 320 ± 390, and 310 ± 660 
Lq yr“^ and with rms residuals of 1400, 3200, 1600, and 
2900 Lq in [/, B, V, and R respectively, where the uncer¬ 
tainties include both statistical uncertainties and systematic 
uncertainties estimated from the rms of light curves of a 
sigma-clipped grid of points within 45 arc-sec of the coordi¬ 
nates of SN 1997bs). The new BST observations reveal that 
the T5551F luminosity at the position of 97bs decreased by 
~ 1.1 ± 0.6 mag between 2001 and 2013 and the TSldlU 
luminosity decreased by ~ 0.7 ± 0.2 mag between 2009 and 
2013. In principle, the BST data can place more stringent 
limits on variability of the source due to the longer baseline 
and better resolution, but in practice the LET variability 
limits may be more reliable and less susceptible to system¬ 
atic errors by virtue of coming from a single instrument and 
using image subtraction. 


2.3 Confusion 

Since we are considering a faint source in a crowded field, we 
must evaluate the likelihood that our detected BST source 
could be an incidental detection of an unrelated source. 
We And that the surface density of all dolphot sources 
within 0”8 of the position of SN 1997bs is 31/arc-sec^, which 
corresponds to a 15% chance of an unrelated source be¬ 
ing detected within 0704 (the distance of the closest source 
from the SN location) by chance. For sources as bright (in 
T8141F) as the detection, the surface density is 7/sq-arcsec, 
which reduces the chance of confusion to 4%. However, aper¬ 
ture photometry seems to indicate that confusion is more 
substantial, with ~ 30% of apertures laid out over a 0”8 grid 
surrounding the source yielding fluxes as bright as our detec¬ 
tion. The large difference between the two methods might be 
due to DOLPHOT limiting detections to only nearly point-llke 
objects or from source deblending. 

Another possibility is that the detection is a surviving 
companion to the progenitor of SN 1997bs. Massive stars 
have a large multiplicity fraction (> 82%) (Chini et al. 
2012). Kochanek (2009) estimates the magnitude distribu¬ 
tion of surviving companions of ccSNe for a set of progenitor 
masses assuming a uniform distribution of mass ratios. Un- 

000-000 


0000 











Imposter SN 1997hs 5 


Table 2. Photometry 


Filter 

Magnitude (± stat ± sys) 

Epoch 

HST WFC3/UVIS F275IF 

> 24.7 

2014-02-08 

HST WFC3/UVIS F336IF 

26.72 ±0.69 ±0.63 

2014-02-08 

HST WFC3/UVIS F4381F 

27.13 ±0.35 ±0.54 

2014-02-08 

HST WFC3/UVIS F5551F 

26.90 ±0.15 ±0.52 

2013-11-28, 2014-02-08 

HST WFC3/UVIS F814IF 

25.49 ±0.12 ±0.11 

2013-11-28, 2014-02-08 

HST WFC3/IR FllOIF 

23.83 ±0.04 ±0.35 

2013-11-28 

HST WFC3/IR F160IF 

22.64 ± 0.04 ± 0.26 

2013-11-29 

SST IRAC 3.6 lira 

> 18.2 (< 15 tjJy) 

July 2013, Feb-Mar-Aug 2014 

SST IRAC 4.5 lira 

> 18.5 (< 7 iiJy) 

July 2013, Feb-Mar-Aug 2014 

SST IRAC 5.8 lira 

> 15.1 (< 109 /iJy) 

May 2004 

SST IRAC 8.0 Atm 

> 12.5 (< 645 /iJy) 

May 2004 

SST IRAC 3.6 lira 

> 18.2 (< 15 tjJy) 

May 2004 

SST IRAC 4.5 lira 

> 16.9 (< 31 iiJj) 

May 2004 

HST ACS/WFC F4351F 

26.78 ±0.17 ±0.04 

2004-12-31 

HST ACS/WFC F5551F 

25.97 ± 0.17 ± 0.03 

2009-12-14 

HST ACS/WFC F8141F 

24.81 ±0.10 ±0.00 

2009-12-14 

HST WFPC2 F555IF 

25.8 ±0.3 “ 

2001-03-04 

HST WFPC2 F814IF 

> 25.0 “ 

2001-02-24, 2001-05-28 


“Photometry taken from Li et al. (2002) 

Magnitudes listed with uncertainties are from DOLPHOT PSF photometry. 
Magnitudes with only upper limits are from aperture photometry. 


fortunately, the mass of the progenitor of SN 1997bs is not 
well constrained, but out of the four progenitor models cal¬ 
culated, the progenitor’s m^sssm = 22.86 is most consistent 
with the 20Mq model. For this progenitor mass and assum¬ 
ing a binary fraction of 80% there is very roughly a 5-20% 
chance that there would be a surviving companion to SN 
1997bs brighter than our detection. However, as discussed 
in Kochanek (2009), most secondaries of exploding stars are 
fainter, blue main sequence stars, but our detection is rela¬ 
tively red (1/ — 7 ~ 1.4) and would correspond to T* ~ 4400 
K if the SED is not strongly influenced by dust. This reduces 
the likelihood that the detection is a surviving companion. 

For the SST data the variability in the 4.5 fim flux 
strongly suggests that the 4.5 fim flux observed in 2004 is as¬ 
sociated with SN 1997bs. The IR flux remaining in 2013/14 
may be only due to diffuse emission coincident with the orig¬ 
inal source. 


2.4 Basic Scalings for Dust 

A common picture of SN impostors is that a surviving star is 
obscured by a dusty shell ejected at the time of the transient. 
For such a model, geometry dictates that at late times the 
evolution of the optical depth of a uniform shell expanding at 
a constant velocity is T{t) oc where t is the time elapsed 
since the ejection of the shell. The observed visible light is 
dominated by scattered photons rather than direct emis¬ 
sion, so inhomogeneities in the shells have less effect than 
naively expected (see Kochanek et al. 2012). In particular, 
inhomogeneities can only grow with time and accelerate the 
optical depth evolution, which would further strengthen our 
arguments. The evolution of the optical depth in turn must 
result in a change in the observed luminosity (in a given 
filter). Consequently, the limits on the observed variability 
of a source surrounded by an expanding shell can be used 
to constrain the current effective absorption optical depth. 


Teff = [Tabs(Tabs + Tsca)]^^^, where Tabs and Tsca are the ab¬ 
sorption and scattering optical depths, by 


1 i f fiLy^ohs \ 
< 2)^ ) ’ 


( 1 ) 


where t is again the elapsed time, Tv.obs is the current ob¬ 
served luminosity, and dLy^obs/dt is the limit on the rate of 
change in the observed F-band luminosity. Similarly, vari¬ 
ability limits constrain the maximum luminosity, 7/*,^, of a 
star surviving within an expanding shell to 


L*,v < 


It/ dLy^obs \ 
2 Tv.eflf \ dt J 


e^v.eft 


( 2 ) 


Late-time data, in addition to constraining the optical depth 
of an expanding shell and the luminosity of a surviving star 
from the observed variability, can also be used to estimate 
the mass of the shell. The mass of the shell, Mej, is related 
to the total optical depth, tv, tot, by 

_ 47rw2tVv,tot(f) 


where Ve is the radial velocity of the shell and Ky is the 
opacity at F-band. The total and effective optical depths 
are related by the albedo w with Teff = (1 — w)^^^Ttot. 

Another case to consider is that the progenitor star is 
obscured by a steady-state dusty wind. We can estimate the 
mass-loss rate needed to obscure the progenitor star with 
such a wind. If we assume that all of the dust forms at the 
dust formation radius, Rf, then the rate of mass-loss is: 


M = 


dTTUwTtif T\/,tot 
Ky 


(4) 


where Rf ~ jTf and Tf ~ 1500 K is the dust formation 
temperature. More precisely the dust begins to form at Rf, 
where the temperature is low enough for grain condensation 
at a rate proportional to density. As discussed in Kochanek 
(2011b), Rf has a complex temperature dependence, but 
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would only range from 2.3 x 10^^ cm for a 7500 

K star to 2.8 x 10^® cm for a 10^'^Lq, 20,000 K star. The 
dust opacity depends on the distribution of grain sizes and 
the dust-to-gas ratio, but will generally be within a factor 
of 2 of 100 cm^/g (of gas) at V-band. If the dust is aspheri- 
cally distributed, the ejected mass inferred by Eqn. 3 for the 
shell case and the rate of mass-loss for the wind case from 
Eqn. 4 could be over- or underestimates depending on the 
specific geometry and orientation of the system. However, 
the temporal scalings for the optical depth encoded by Eqn. 
1 and 2 should remain reasonable approximations even for 
complex geometries (Kochanek et al. 2012). 

We will utilize these relations in §3.1 and §3.2 to help 
determine whether a surviving star to SN 1997bs could be 
obscured by an expanding shell or steady-state wind. 


2.5 DUSTY 

We model the SED of the source using dusty (Ivezic & 
Elitzur 1997; Ivezic et al. 1999; Elitzur & Ivezic 2001), a 
code for solving radiative transfer through a spherically sym¬ 
metric dusty medium. We use stellar atmospheric models 
from Castelli & Kurucz (2004) for stars of various tempera¬ 
tures and solar composition. We find best-fitting models us¬ 
ing a Markov Chain Monte Carlo (MCMC) wrapper around 
DUSTY for both silicate and for graphitic dust from Draine 
& Lee (1984) using a standard MRN grain size distribution 
{dn/da oc ® with 0.005 /rm < a < 0.25 /rm; Mathis et al. 
1977). In Appendix A we show that our conclusions are not 
dependent on this choice for the grain size distribution. Since 
silicate dust is the type expected to form around more mas¬ 
sive stars (like the progenitor of SN 1997bs was believed to 
be) and the conclusions of the paper are robust to the cho¬ 
sen dust type, we will only present the results for silicate 
dust. For the MRN grain size distribution of silicate dust 
wv — 0.86, making rv.eff — 0.37Tvyot. The models assume 
that a shell with a thickness (Rout/Rin) of 2.0 expanding 
with a velocity of Ve = 765 kms“^ (± a factor of 2; Smith 
et al. 2011) was ejected during the 1997 event. We actively 
expand Rin at this rate since the effects of dust are domi¬ 
nated by Rin because optical depth drops as R“^. We also 
produce models for a dusty wind by setting Rin equal to the 
radius corresponding to a typical dust formation tempera¬ 
ture of 1500 K (see, e.g., Kochanek 2014b). The used in 
the MCMC is calculated with the logarithmic (linear) dif¬ 
ferences between modeled and observed fluxes (or limits) in 
each filter treated as a detection (or limit) and, for the ex¬ 
panding shell models, between modeled and observed Ve (to 
place the dust at the Rin corresponding to the given dust 
temperature). When considering all photometric constraints 
as only upper limits, we calculate the Ax^ compared to hav¬ 
ing no star as functions of R* and T*. When folding in the 
light curves we add x^ contributions for each observation, 
n, in each Hlter, /, found by 


2 / Lohs.f.n Lmod,f,n \ 

Xf,u 

(5) 

V J 


where the model luminosity 


Rmod,/,n = L,je 

(6) 


is controlled by the evolution of the optical depth 



Figure 3. The SED of SN 1997bs with the best-fitting model for 
a silicate shell when treating our F555W and F814iy photometry 
as detections and the other bands as upper limits. The latest HST 
and SST detections and limits for a surviving star are shown as 
the open pentagons with the current best-fitting observed model 
spectrum shown as the thick black line. This best-fitting model 
has Tt = 18,900 K, ry tot = 5.34, and L, = 10'*'®'^ Lq. We also 
show this best-fitting model evolved forwards in time to 2017 in 
green and backwards in time to 2009 in blue with r oc but 
with the same T* and L*. The tension between the 2009 HST pho¬ 
tometry, shown as the open blue triangles, with the best-fitting 
shell model evolved back to 2009 illustrates how models in which 
SN 1997bs is obscured by an expanding dusty shell cannot be rec¬ 
onciled with the late-time optical evolution. Although the present- 
day luminosity constraints allow for SED solutions with higher IR 
fluxes these solutions require higher optical depths, which, when 
evolved backwards in time, are even more incompatible with the 
2009 HST photometry. Eor comparison, we show the HST pro¬ 
genitor detection as the solid black square and an unobscured 
model spectrum with T* = 18,900 K and L* = 10® '*^ Lq as the 
dashed magenta line. 







(7) 


In this model, t^sj is the optical depth at the time of the 
latest observations T = 16.5 yr after the transient and t„ is 
the time elapsed since the transient for observation n. 


3 A STAR OBSCURED BY DUST? 

We consider whether the progenitor of SN 1997bs could have 
survived obscured by either an expanding dusty shell or by 
a steady-state wind. Our candidate late-time detection, like 
those in earlier studies, is marginal and confusion is a non¬ 
trivial issue. Therefore, we will consider a range of inter¬ 
pretations for which observations constitute detections and 
which are only upper limits. We consider the case of F555W 
and R814TU as detections since detections in these filters 
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have been previously reported (Li et al. 2002; Van Dyk & 
Matheson 2012). We also consider the case of only F8141V as 
a detection, as this is the filter in which a detection is most 
convincing visually. We consider the case of F81AW and 
4.5 fim as detections because there was a decline in 4.5 /im 
flux between 2004 and 2014 and this represents a scenario 
in which there is high obscuration. Finally, we consider the 
case where all of the photometry is treated as upper limits, 
since the candidate detections are marginal and could be 
due to confusion. 


3.1 Obscuration by an Expanding Shell 

First, we model our photometry as though a surviving star 
obscured by an expanding shell is recovered with detections 
in EdhSlF and ESldlF but only upper limits in the other 
Hlters. We present, as an example, the best-htting SED for 
this case in Fig. 3. This best-fitting SED is almost a factor of 
4 fainter than the luminosity of a progenitor with the same 
T* and relies on a large amount of obscuration {tv, tot = 5.3) 
to account for the low optical flux currently observed. 

The complete results of our MCMC modeling are shown 
in Fig. 4. When considering only the latest photometric con¬ 
straints for each hlter (the cases displayed on the left side 
of Fig. 4), it is the possible that the star survived with 
unchanged luminosity but is cloaked beneath a significant 
amount of extinction. In the cases where either the EhhSlT 
and F814W photometry or only the ESldVF photometry is 
treated as a detection, the star can be cool with relatively 
low luminosity and low obscuration or hot with higher lu¬ 
minosity and higher obscuration. The E814VF-only case al¬ 
lows a more luminous surviving star than the ^55514^ and 
E8144F case because the models no longer have to ht the 
shallow slope between the F’555W and ^81414^ magnitudes, 
enabling scenarios with higher optical depths to achieve 
good fits. In the F'8144F and 4.5 /rm case, a cooler star with 
lower luminosity and lower obscuration is not allowed by the 
mid-IR detection, which requires a significant luminosity to 
be reprocessed by dust. If all of the photometry is taken 
as only upper limits, the MCMC modeling would be poorly 
constrained, so we instead show the maximum allowed lu¬ 
minosity of a surviving star for optical depths of tv, tot = 
0, 1, 3, and 10. The limits-only case echoes the results of 
the other cases - the latest photometric constraints, taken 
alone, allow for a luminous, but heavily obscured, surviving 
star. 

These limits are relatively robust to variations in the ve¬ 
locity of the expanding shell. Doubling the assumed velocity 
of the expanding shell only significantly increases the maxi¬ 
mum luminosity of a surviving star with T* > 25, 000 K and 
Tv, tot ~ 10 and then only by ~ 0.3 dex. Conversely, halving 
the expansion velocity decreases the maximum luminosity 
in the same parameter range by a similar amount. 

That solutions consistent with the photometric data ex¬ 
ist does not mean they are physical because they may im¬ 
ply unphysical ejecta masses or unobserved variability. First, 
we consider the ejecta mass needed to obscure a surviving 
star with an unchanged intrinsic luminosity and tempera¬ 
ture. Fig. 4 shows that this would require the current optical 


depth to be Tv,tot ~ 10. Using Eqn. 3, we estimate 



Such an ejected mass is feasible, but, as we shall see, this sce¬ 
nario is inconsistent with the variability constraints dictated 
by geometric expansion if any of the late-time photometric 
detections are powered by the luminosity from a surviving 
star. The kinetic energy of the ejecta, in 

this scenario would be 



We compare this to the radiated energy, Erad = U.sC^peak, 
of the transient estimated by Smith et al. (2011) of 


where Lpeak is the peak luminosity of the outburst, ti .5 is the 
time for the transient to fade by 1.5 mag from its peak, and 
^ ~ 1 is a dimensionless factor that depends on the shape of 
the light curve. The high ratio of radiated to kinetic energy 
in this scenario would likely require a radiative, rather than 
explosive, mechanism. We do caution that these estimates of 
Mej and Ekin assume homogeneous, spherically distributed 
ejecta. 

Fig. 3 also illustrates how the expanding shell model 
that best hts the latest data would necessarily have been 
much fainter in the optical in 2009, in gross disagreement 
with the HST photometry from 2009. A model in which 
Tv,tot = 10 would be in even more severe conflict with the 
photometric evolution. More quantitatively, using Eqn. 1 
with the 3cr upper limit on the variability from the LET 
data constrains the current maximum optical depth of an 
expanding shell to be 


TV,eS < 11 


/ t \ / 1140 Lq \ f dLv,ohs/dt \ 
1^16.5 yr J ^YjOhs J 1500 Lq yr“i J 


Similarly, we can use the variability limits from LET and 
Eqn. 2 to constrain the maximum luminosity of a star sur¬ 
viving within an expanding shell as 


L*.v < 3.3xl0'‘ 


t 


dLv, ohs /dt 

16.5 yr / V 1500 Lq yr~ 


0'^ V.etf — 1 

TV,ea 


-Lq.{12) 


These limits from the LET are less constraining than the 
decline in flux seen in the HST observations. In Fig. 4, the 
panels on the right show the constraining effect of using the 
both the LET and HST light curves from 2008-14 with Eqn. 
5 and 6. The effects of the luminosity constraints should 
be relatively robust to departures from our assumptions of 
homogeneous, spherically-distributed dust. The light curves 
restrict solutions to have low optical depths (rv,tot < 1 ). 
Although the 7-1-4.5 fim+LC case still appears to favor a 
luminous, obscured survivor, the Xmin for this case is much 
higher than that for the other cases, indicating that it is 
strongly disfavored. Moreover, the putative HST source at 
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Figure 4. MCMC results for the luminosity and temperature of a surviving star obscured by an expanding shell ejected at the time of 
the SN 1997bs transient when imposing different possible photometric constraints: treating the HST WFC3/UVIS F555VF and FSIAW 
photometry as detections and the other filters as only upper limits (top row), treating the HST WFC3/UVIS F814VF photometry as 
a detection and the other filters as upper limits (2nd row from the top), treating the HST WFC3/UVIS F814VF and SST 4.5 ^m 
photometry as detections and the other filters as upper limits (2nd row from the bottom), and treating all photometry as upper limits 
(bottom row). The panels on the right also fold in the constraints from the LBT and HST light curves. The shaded bands in the panels 
on the bottom row show the luminosities within the 90-99.99% confidence intervals (6.25 < < 21.1 for three parameters — L*, 

T*, and r) relative to no surviving star) for r\/tot = 0^ 1? 3, and 10. For comparison, the solid black line and gray band indicates the 
progenitor luminosity as constrained by the pre-explosion measurement of mpeoe = 22.86 and its la uncertainty (0.16 mag). The lowest 
value for the 10,000 MCMC steps accepted is given for each panel. The panels on the left show that the latest photometric constraints 
taken alone do not rule out the star surviving with its pre-eruption luminosity if obscured by a shell with ry^tot > 5. However, the panels 
on the right that fold in the light curve evolution only allow a luminous surviving star in the ‘/-|-4.5^m-|-LC’ and ‘Limits-|-LC’ cases. 
The former is heavily disfavored by the relatively high X^in the latter would require that all late-time optical detections are of an 

unrelated source. 


0000 


000 , 


000-000 


















Imposter SN 1997bs 9 


the coordinates of 97bs seems to have faded between 1998 
and 2013, which further disfavors any expanding shell sce¬ 
nario in which the optical flux at the position of 97bs is from 
a surviving star. The remaining expanding shell scenario is 
that the star is still heavily obscured (rv.tot > 10) and the 
optical emission reported here and by Li et al. (2002) and 
Van Dyk & Matheson (2012) is not from a surviving star. 

Such a scenario is difficult to reconcile with the hrst 
year of the transient light curve. The light curve presented 
by Van Dyk et al. (2000) shows that SN 1997bs rapidly 
evolved to the red, with the redward color evolution peaking 
at V —7 ~ 3 by eight months after discovery, which given the 
constraints on the ejecta velocity and the progenitor V-band 
detection is roughly the time at which the ejecta should have 
reached the dust formation radius. As discussed in Kochanek 
et al. (2012), this peak in the color evolution corresponds 
to a peak in the optical depth at tv — 10. A shell with 
Tv ~ 10 at eight months would have evolved to tv ~ 0.02 
by our latest HST observations. Such a low optical depth 
would not allow a luminous survivor. It may be possible that 
slower-moving ejecta reached the dust formation radius and 
formed additional dust during gaps in the coverage of the 
light curve later than nine months post-discovery, but all 
subsequent observations have constrained SN 1997bs to be 
much less red than it was at eight months. 


3.2 Obscuration by a Wind 

Rather than being obscured by a dusty shell ejected during 
the outburst in 1997, the surviving star could be obscured by 
a steady-state dusty wind that began following the transient. 
For our standard models we adopt an inner (dust formation) 
temperature of Tf = 1500 K. We present, as an example, the 
best-fitting SED for a wind when treating our 7^555IV and 
7^814tV photometry as detections and the other bands as 
upper limits in Fig. 5. The striking fact is that the luminosity 
of the best-fitting SED for our candidate detection is fainter 
than the progenitor by nearly a factor of 20 when modeled 
as being obscured by a dusty wind. 

The complete MCMC results are shown in Fig. 6. The 
HST FllOW and F160W photometric limits strongly con¬ 
strain the amount of hot dust around the source. As a result, 
at least for the cases where only optical filters are treated 
as detections, a surviving star is constrained to be cool with 
low to moderate obscuration {Tv,tot < 7) and a significantly 
lower luminosity than the progenitor. While the case where 
the 7^81411^ and 4.5 pm photometry are treated as detections 
appears to favor a heavily obscured source with only slightly 
lower luminosity than the progenitor, all of the models in 
this case have high values because of the great difhculty 
in accounting for the 4.5 pm flux with T{ = 1500 K without 
violating the HST near-IR limits. When we model all of the 
photometry as upper limits (see the lower-left panel in Fig. 
6 ) we see that tv, tot > 10 is needed to allow the star to sur¬ 
vive with an undiminished luminosity behind a dusty wind. 
As illustrated by the high values of the case where the 
7^814IF and 4.5 pm photometry are treated as detections, 
such a high optical depth is only possible if all the optical 
constraints are only upper limits. 

These limits are relatively robust to variations in the as¬ 
sumed dust condensation temperature in the wind. Models 



Figure 5. The SED of SN 1997bs with the best-fit model for a 
silicate wind with an inner edge temperature of Tf = 1500 K when 
treating our F555W and i^8141T photometry as detections and 
the other bands as upper limits. The HST and SST detections 
and limits for a surviving star are shown as the black open pen¬ 
tagons with the best-fit observed model spectrum shown as the 
thick black line, the intrinsic input SED is the thin blue line, and 
the contributions of the attenuated input radiation, scattered ra¬ 
diation, and dust emission are shown as thin red solid, dotted, and 
dashed lines respectively. This particular model has T* = 5252 
Tv, tot = 1.95, and T* = 3197 Lq. While the mid-IR limits 
are relatively weak, the near-IR constraints on hot dust prevent 
solutions with much higher luminosities and optical depths. For 
comparison, we show the HST progenitor detection as the solid 
black square and a model spectrum also with T* = 5252 K and 
L* = 10^ ®'^ Lq as the dashed magenta line. 


of a steady-state wind with cooler dust condensation tem¬ 
peratures result in slightly weaker luminosity limits for a 
hot, heavily obscured surviving star. Even if we drop the 
condensation temperature to If = 1000 K (750 K) the lumi¬ 
nosity limit for a T* > 7500 K surviving star with tv, tot = 10 
only increases by ~ 0.3 dex (~ 0.5 dex). We note, however, 
that setting Tf = 1000 K significantly improves the fit of the 
model (Xmin = 0.9) in the case where 7^814W and 4.5 /rm are 
treated as detections by easing the tension with the near-IR 
flux limits. 

We can use Eqn. 4 to estimate the mass-loss rate re¬ 
quired to achieve tv, tot ~ 10. Using the observed ejecta ve¬ 
locity of SN 1997bs as Uw, the mass-loss rate would corre¬ 
spond to 3x lO”'* MQyr“^ (4x 10“® Mqjx~^) for a 10“^'^ Lq 
(10® '* Lq), 7500 K (20,000 K) star. Such high mass-loss rates 
are typical of active LBVs (Puls et al. 2008), although a still 
higher mass-loss rate of > 10“^'® Mq yr“* is needed in or¬ 
der for dust formation to occur around hot (T > 10,000 K) 
stars (Kochanek 2011b, 2014b). Additionally, at these high 
mass-loss rates, the wind itself becomes optically thick in¬ 
dependent of any dust leading to a pseudophotosphere with 
an apparent temperature of T ~ 7000 K (Davidson 1987). 
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Figure 6. MCMC results for the luminosity and temperature of a surviving star obscured by a steady-state wind when imposing different 
possible photometric constraints: treating the HST WFC3/UVIS F555W and F81AW photometry as detections and the other filters as 
only upper limits (top left), treating the HST WFC3/UVIS i*’814VF photometry as a detection and the other filters as upper limits (top 
right), treating the HST WFC3/UVIS and SST 4.5 /rm photometry as detections and the other filters as upper limits (bottom 

left), and treating all photometry as upper limits (bottom right). The shaded bands in the ‘Limits Only’ panel show the luminosities 
within the 90-99.99% confidence intervals (6.25 < < 21.1 for three parameters — L*, T*, and r - relative to no surviving star) 

for TV,tot = 0, 1, 3, and 10. For comparison, the solid black line and gray band indicates the progenitor luminosity as constrained by 
the pre-explosion measurement of rriFeoe = 22.86 and its Itr uncertainty (0.16 mag). The lowest value for the 10,000 MCMC steps 
accepted is given for each panel. Well-fitted wind models are possible when treating the 1^55514^ and F814W photometry as detections, 
but these models require a surviving star to be much fainter than the progenitor. Though the panel where both the F81AW and the 
4.5 /rm photometry are treated as detections appears to allow a luminous survivor, these models have a high x^ because of the great 
difficulty in accounting for the 4.5 fim flux with Tf = 1500 K dust without violating the HST near-IR limits. It is unlikely that a surviving 
star is obscured by a steady-state wind unless the star is significantly fainter than the progenitor. 


A mass loss rate of M > 10 Mq yr ^ would cor¬ 
respond to TV,tot > 140 (12) for r* = 7500 K (20,000 K). 
Essentially this means that for a hot surviving star there 
must be either no dust in the wind or a very thick dusty 
wind {tv, tot > 12). However, Tv,tot >> 10 is inconsistent 
with any optical detection of a source. 


4 DISCUSSION 

The new HST images we present clearly lack a source with 
optical flux comparable to the progenitor of SN 1997bs. If 
SN 1997bs was a non-terminal event the diminished optical 
flux could be due to increased obscuration of the star, an 
increase in the star’s temperature, or an intrinsic decline 
in the luminosity. The alternative is that SN 1997bs was a 

000 , 


genuine SN. We evaluate these different scenarios and the 
limits placed on them by the data. 


4.1 Obscuration 

The late-time photometric evolution of SN 1997bs poses a 
challenge to the prevailing notion of it being an SN im¬ 
postor in which a non-terminal eruption ejected a signifi¬ 
cant amount of mass that then obscured the surviving star. 
The observed limits on the variability at the position of SN 
1997bs from HST and LET rule out any expanding dusty 
shell scenario in which the optical emission represents a ‘re¬ 
covery’ of the luminous progenitor (see the right-hand panels 
of Fig. 4). The data cannot, however, rule out the possibility 
that a surviving star is still cloaked by a dusty shell with 
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Tv,tot >10, but this scenario is difficult to reconcile with the 
earlier color evolution of the transient. 

The new data are also difficult to reconcile with the al¬ 
ternative, a steady dusty wind, which Kochanek et al. (2012) 
found was consistent with earlier archival data (see Fig. 6). 
The high mass loss rate needed for self-shielding so that dust 
can form in the wind of a hot (T* > 10, 000 K) star would 
result in a large tv, tot that would not allow the candidate 
optical detections reported here and by Li et al. (2002) and 
Van Dyk & Matheson (2012). 


4.2 Bolometric Correction 

If we assume that there is no dust obscuration, the maxi¬ 
mum luminosity of the surviving star is ~ lO"*'^ Lq (see the 
Tv,tot = 0 line in the ‘Limits Only’ panel of Fig. 4 or 6), 
which is significantly fainter than the traditional range of 
luminosities for LB Vs. It is unlikely that the low luminosity 
could solely be the result of a large bolometric correction. 
Allowing a surviving star to have an effective temperature 
of up to 40,000 K, the maximum luminosity of a surviving 
star is still at least a factor of three fainter than the mini¬ 
mum luminosity of the progenitor. For temperatures higher 
than 40,000 K, the maximum luminosity of the star can 
be constrained by the maximum possible H/3 flux allowed 
by our TSSSIT photometry, since a very luminous hot star 
would produce a significant flux of ionizing photons and the 
ejecta would provide an absorbing medium. If we assume 
that the star emits a blackbody spectrum and is surrounded 
by enough gas to absorb the photons with 11% of recombina¬ 
tions going through the H/3 branch (based on Case B recom¬ 
bination coefficients for T* = 20, 000 K from Draine (2011)) 
and set the H/3 flux equal to our flux limit for f’SSSIF (we 
do not use Ha because its wavelength places it in the tail 
of the TSldlF filter response curve), we find that the maxi¬ 
mum luminosity of a T* = 40, 000 K star is ~ 10®'° Lq, with 
lower maximum luminosities for increasing T,. Accordingly, 
we cannot completely rule out the possibility that a pro¬ 
genitor initially with 4400 < T* < 13, 000 K survived the 
episode as a T, > 40,000 K Wolf-Rayet star. However, we 
are unaware of any theoretical studies suggesting that a sin¬ 
gle short duration event could affect such a dramatic change. 


4.3 A Tuckered-Out Star 

Another possibility that would allow SN 1997bs to be a non¬ 
terminal event is that the progenitor survived the event with 
a diminished luminosity - the ‘tuckered-out’ star hypothesis 
(Smith et al. 2011). If a surviving star is unobscured and its 
temperature is unchanged from pre-emption, it must be at 
least a factor of ~ 30 fainter than the progenitor. No known 
mechanism could so drastically diminish the intrinsic lumi¬ 
nosity of the star. As mentioned in the previous section, even 
if we allow for a significant increase in the temperature of the 
star, a factor of four decrease in the luminosity is required 
by the latest observations. A conspiracy of factors could ex¬ 
plain the dramatic drop in apparent luminosity. Though an 
expanding dusty shell is ruled out by limits on the observed 
variability, a combination of a steady-state dusty wind and 
a significant increase in T* could allow for a surviving star 
with a luminosity that has declined by as little as ~ 25%. 


However, this collusion is disfavored by the candidate late¬ 
time detections of an optical source. 

The ‘buildup’ or ‘recovery’ time-scale for the radi¬ 
ated energy budget given by Smith et al. (2011) is trad = 
ifrad/L. = ti.sC-t'peak/L*, where Flrad is the energy radiated 
during the outburst, Lpeak is the peak luminosity of the 
outburst, L* is the quiescent pre-outburst luminosity of the 
progenitor, and C is a dimensionless factor. For SN 1997bs 
this is 


which means that a large fraction of this recovery time-scale 
has already elapsed with no re-brightening of a surviving 
star. 

This time-scale may not be the most relevant time-scale 
for the stellar luminosity, since the envelope would likely 
return to thermal equilibrium primarily through Kelvin- 
Helmholtz contraction rather than by energy radiated from 
the core. In order to decrease the bolometric luminosity for 
such an extended period of time, either some of the nuclear 
energy generated must be diverted into gravitational poten¬ 
tial energy by inflating the stellar envelope or the nuclear 
luminosity must have decreased. 

If the missing luminosity is being used to alter the struc¬ 
ture of the outer layers of the star, the time-scale for it to 
lift mass AM against the gravitational potential of the star 
is 


tinfl ~ 120 


M* 


V M© J y 20 Mq 
L* 

10® Lq 


10,000 K 

- 3/2 


Offi 

T 


yr, 


(14) 


where M*, T*, and L* are the mass, temperature, and lu¬ 
minosity of the progenitor, and /, is the fraction by which 
the radiated luminosity has decreased. This could account 
for the missing energy, but it is unclear if there is a the¬ 
oretical mechanism that could drive this. Moreover, LoR is 
not so long that we would expect no external changes on 
time-scales of a decade. Another problem is that the most 
natural post-eruption state would be a star with an over- 
expanded envelope rather than the reverse. The envelope 
would then shrink on a thermal time-scale, rapidly at first 
but then slowing, making the star overluminous, not sublu- 
minous. The overexpansion is a natural consequence of any 
transient mechanism which has no ‘knowledge’ of the es¬ 
cape speed and has mainly been discussed in the context of 
shock-heating non-degenerate companions of Type la SNe 
(Pan et al. 2013; Shappee et al. 2013). 

The final alternative is for the nuclear luminosity to 
have decreased. The core luminosity should be unaltered as 
it can only change over a (core) thermal time-scale, but it 
might be possible for the luminosity from shell-burning to 
decrease suddenly in certain (fine-tuned) situations. Perhaps 
the progenitor had previously experienced enough mass loss 
for its hydrogen-burning shell to be close to its surface and 
the star experienced a final shell flash, analogous to the final 
shell flash of an asymptotic giant branch star. If the progen¬ 
itor was very massive star (Mi > 80 Mq), a thermonuclear 
outburst from pulsational pair instability might be another 
possibility (Woosley & Heger 2015). 
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4.4 SN and Confusion 

Alternatively, SN 1997bs could have been a terminal event. 
If the star went supernova, our detection of an F9>l'iW = 
25.49 ± 0.12 source 07016 from the position of SN 1997bs 
must be due to confusion, a surviving companion, or the SN 
ejecta. As we describe in §2.3, the likelihood of our detection 
being due to confusion is between 4% and 30%. The likeli¬ 
hood of the source being a surviving binary companion is 
comparable, at 5-20%, although we would expect a surviv¬ 
ing companion to be bluer. The slight decrease in luminosity 
between the 2001 i/'S'r/WFPC2, 2009 HST/ACS, and 2013 
HST/WFC3 F555W and F814IT images reduces the likeli¬ 
hood that the detection is due to confusion, although there 
could be systematic issues in comparing the crowded field 
photometry based on the three different instruments and the 
LET variability limits are consistent with zero. 

The 4.5 /rm flux measured at the location of SN 1997bs 
dropped from a significant 31 ± 4/iJy in 2004 to 7 ± 3/iJy 
in 2014. This variability indicates that the 4.5 /rm flux ob¬ 
served in 2004 originated from SN 1997bs and is not confu¬ 
sion. The decrease in the 4.5 /rm flux while the 3.6 /rm flux 
stayed constant is difficult to explain if the 3.6 /im flux is 
from SN 1997bs. More likely, the 3.6 /im flux is unrelated 
to SN 1997bs (SN 1997bs is close to a spiral arm) and the 
decreasing 4.5 /rm emission was from warm (~ 1000 K) dust 
that has since cooled. 

It is possible that we are detecting residual flux from 
the shock interaction of SN 1997bs with its circumstellar 
medium (CSM). If we parametrize the fraction of the max¬ 
imum possible shock luminosity radiated in our observed 
hlters as /, then the observed shock luminosity is given by 


Ls,obs ^ 2.3 X 10^ 


M 


IO-^Mq 


;)( 


Ve 


765 kms-i 
^ / 


(lool^)’ 


where M is the pre-eruption mass loss rate, Ve is the velocity 
of the SN ejecta, and Vw is the pre-eruption wind velocity. 
Although we do not have constraints on M or since little 
is known about the progenitor, reasonable values of and 
/, together with a relatively (but not outrageously) high M 
during the decades proceeding the transient, are sufficient to 
possibly account for our detection (~ 1700 Lq in FSl'iW). 

It is worth reconsidering why SN 1997bs was originally 
designated as a likely SN impostor: a low peak luminosity 
and a possible flattening of the late-time (> 250 days post¬ 
explosion) light curve at ~ 0.5 mag fainter than the progen¬ 
itor (see Van Dyk et al. 1999). Though the source continued 
to fade well beyond this level. Van Dyk & Matheson (2012) 
re-asserted the case for a surviving star based on the 2009 
HST observations in which the star was ‘recovered’ with 
a brightness consistent with measurements made in 2001 
(see Van Dyk & Matheson 2012). However, the source at 
the coordinates of SN 1997bs, rather than re-brightening as 
required if the optical emission was from a surviving star 
partially obscured by an expanding dusty shell, has only 
continued to fade and is now ~ 4 mag fainter than the pro¬ 
genitor. For the most likely physical conditions, the near-IR 
limits rule out the detection of optical emission from a star 
surviving in a dusty wind. Rather than providing evidence 


that 97bs was a non-terminal eruption, the late-time obser¬ 
vations are difficult to reconcile with this prevailing notion. 

SN 1997bs can be explained as a low-energy SN. Over 
the last two decades, a class of faint Type IIP SNe has 
emerged with peak magnitudes as faint as My = —13.76 
(SN 1999br; Pastorello et al. 2004), less than a magnitude 
brighter than 97bs {My ~ —12.9). These faint SNe have 
low ejected ®®Ni masses (~ 10“^M©) and require low ex¬ 
plosion energies [L < 10®^ erg). The low ®®Ni mass is likely 
due to fall-back of material on to the collapsed remnant 
(see e.g., Woosley & Weaver 1995). Although the late-time 
light curve of SN 1997bs might have significant contributions 
from a CSM-SN shock interaction, we can estimate the max¬ 
imum ®®Ni mass ejected by scaling the V-band luminosity 
to that of the well-studied faint IIP SN 1997D at the same 
epoch (270 days post-explosion). Nominally this would sug¬ 
gest that that the maximum ®®Ni mass for 97bs is 20 times 
smaller than the 0.002 Mq estimated for SN 1997D by Tu- 
ratto et al. (1998). However, 97bs was much redder than 
97D by this point {V — I ~ 4.5 versus ~ 1.2 mag) and was 
likely significantly obscured. If we posit that extinction is 
needed to match the colors, then 97bs had E{B — V) ~ 2.3 
mag at this point, which would increase the maximum ®®Ni 
mass by nearly a factor of 1000. Even if the shock luminos¬ 
ity contributed a large fraction of the luminosity, 97bs could 
still potentially have more ejected ®®Ni than a significant 
fraction of faint HP’s. 

Thus, SN 1997bs could be a Hn analog of a faint IIP. 
Progenitor detections and hydrodynamical modeling of light 
curves point to faint HP SNe progenitor masses of 10-15 Mq 
(Spiro et al. 2014). Meanwhile, recent work in SN theory by 
Ugliano et al. (2012) and Pejcha & Thompson (2015) sug¬ 
gests that stars with masses between 14 and 16 Mq, some 
between 17 and 22 Mq, and stars between 23 and 26 Mq 
are the most difficult to explode. Faint IIP SNe might arise 
from the lower mass (14-16 Mq) stars that have density 
structures resulting in low-energy explosions while 97bs and 
similar ‘impostors’ could arise from the higher mass (23-26 
Mq) stars. These higher mass stars have the potential to 
have the pre-explosion mass loss required for a transient to 
have the narrow hydrogen emission lines that are character¬ 
istic of SN Iln (and SN impostor) spectra. 

The event could even be a collapse-induced thermonu¬ 
clear SN (Kushnir & Katz 2014), which would dump enough 
energy in the stellar envelope quickly enough to drive a shock 
in a manner similar to the piston simulations of Dessart et al. 
(2010). This would be a weak, terminal event of a star in one 
of the harder-to-explode mass ranges. 

All of these low-energy explosions are also much more 
favorable for dust formation because condensation will occur 
at much higher densities (see, e.g., Kochanek 2014a). This 
would provide a natural explanation of the increasingly red 
colors observed by Van Dyk et al. (2000), since dust for¬ 
mation fairly naturally correlates with dropping luminosity. 
The lack of a surviving star then avoids the problems with 
dusty models of a surviving star. 


4.5 Summary 

Our late-time observations of SN 1997bs seem most easily 
understood if SN 1997bs was, in fact, the hrst SN found by 
LOSS. It is difficult to explain the data as dust obscuration 
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by a wind and it is only possible to explain by an expanding 
shell if the star is still completely obscured in the optical 
and the putative recovery of a surviving star is actually un¬ 
related emission. Similarly, invoking bolometric corrections 
or simply making the star less luminous seems difficult to 
support. While the observations do not rule out a surviv¬ 
ing star, they do undermine the evidence in favor of this 
scenario. Occam’s razor, if nothing else, argues for simply 
making it an SN rather invoking a more baroque model. 

Nature does not, of course, have to respect Occam’s ra¬ 
zor. Further observations are necessary to confirm the fate of 
SN 1997bs. Unfortunately it is impossible to resolve any IR 
emission from SN 1997bs with the diffuse emission from the 
nearby spiral arm with SST. The improved resolution and 
sensitivity of the James Webb Space Telescope are needed 
to reduce the mid-IR limits enough to completely rule out a 
cool star behind a thick, dusty wind or shell. 

As argued by Kochanek et al. (2012), the current class 
of SN impostors is likely comprised of multiple phenomena. 
We find that SN 1997bs, an archetypal SN impostor, was 
likely a genuine SN, but there are also several SN impostors 
that have clearly undergone multiple non-terminal outbursts 
(e.g., SN 2000ch, 2002kg, 2009ip). Moreover, in newly ob¬ 
tained HST images, we confirm a surviving star to the SN 
impostor 1954J, as proposed by Smith et al. (2001) and Van 
Dyk et al. (2005) (Adams & Kochanek, in preparation). At 
least a significant fraction of SN impostor progenitors are 
relatively low mass (15 Mq < M < 25Mq) stars (Smith 
et al. 2011; Kochanek et al. 2012). These masses are lower 
than the mass range of ‘classical’ LBVs and should not be 
close to their Eddington limits and any of the radiative in¬ 
stabilities of classical LBVs. Recurrent impostors have the 
highest progenitor masses (see Table Al). Perhaps some of 
the lower mass impostors are faint Iln SNe and only the 
impostors arising from more massive stars are non-terminal 
eruptions of LBVs. 
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Table Al. SN Impostor Masses 


Object 

M,[Mq] 


Single Episode 


SN 2008S 

~ 10“ 


2008-OT 

12 - 15“ 


SN 1997bs 

12 - 24“, 

> 20'’ 

SN 2003hm 

14- 17“, 

> 20'’ 

SN 1954J 

14-21“, 

> 20'’ 

2009-OT 

~ 25'’ 


HD 5980 

58 - 79“* 


SN 1961V 

> SO'' 


Reccurent 


P Cygni 

20 - 40“* 


SN 2002kg 

> 15“, > 

40'’ 

SN 2000ch 

20 - 80“, 

> 407 

SN 2009ip 

50 - 80'’ 


T) Car 

~ 1609 



‘Current mass 


“Kochanek et al. (2012) 

*’Van Dyk & Matheson (2012) 

“Foellmi et al. (2008) 

‘^Kochanek et al. (2011) 

“Lamers et al. (1983) 

7Wagner et al. (2004) 

^Davidson & Humphreys (1997) 

University of Notre Dame, the University of Minnesota, and 

the University of Virginia. 


APPENDIX A: DUST SIZE DISTRIBUTION 

Although the MRN dust size distribution we assumed for 
the analysis presented in this paper is a reasonably accurate 
prescription for the size distribution of interstellar dust, dust 
formed around a particular star will grow to a particular size 
that is dependent on the density and radiation field of the 
circumstellar environment (see, e.g., Kochanek 2014b, for 
a discussion). Unfortunately, we are unable to determine a 
most likely dust size, since the size, an input to our models, 
is dependent on the output (i.e., r, T*, and L*) from the 
models. 

Instead, we evaluated our models for a number of dust 
sizes ranging from 0.005 fim to 0.25 pm and verified that 
our conclusions are robust to the chosen dust size. For an 
example, in Fig. Al we present the results of the silicate 
wind model for a range of dust sizes and effective opti¬ 
cal depths. For silicate dust grain sizes, a, of 0.005 pm, 
0.035 pm, 0.25 pm, and an MRN size distribution the albedo 
is wv — 0, 0.10, 0.43, and 0.39, corresponding to Tv,eB — 
1.00, 0.82, 0.32, and 0.37, respectively. For a given effective 
optical depth the luminosity limits on a surviving star are 
largely unchanged. The most significant difference is that 
the of fho best-fitting model for a silicate wind with our 
ESld and 4.5 pm photometry treated as detections for a 
dust size of 0.035 pm is only 1.3 (compared to 3.9 for an 
MRN distribution of dust), slightly reducing the confidence 
with which we can rule out a luminous star surviving behind 
a dusty wind. 

The impacts of the dust size on the results of the shell 
model are even more limited. For a given effective optical 
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Figure Al. The luminosities within the 90-99.99% confidence intervals (6.25 < < 21.1 for three parameters - L*, T*, and r - 

relative to no surviving star) for rytot = 0) 1) 3, and 10 with a steady-state silicate wind for different dust grain sizes, a, as compared 
to the results for an MRN dust size distribution (bottom-right panel) when treating all photometry as upper limits. The solid black 
line and gray band indicates the progenitor luminosity as constrained by the pre-explosion measurement of m^goe = 22.86 and its Itr 
uncertainty (0.16 mag). For most temperatures and effective optical depths, using a different (reasonable) dust grain size distribution 
would only shift the luminosity limit by ~ 0.1 dex. 


depth, the only significant change for different dust sizes is 
that the luminosity limit for a cool star with heavy obscn- 
ration is increased. The key point, however, is that the light 
curve constraints continue to rule out the optical emission 
being a detection of a luminous surviving star dimmed by 
an expanding shell. 

REFERENCES 

Alard C., 2000, A&AS, 144, 363 
Alard C., Lupton R. H., 1998, ApJ, 503, 325 
Berlin E., 2006, in Astronomical Society of the Pacific Con¬ 
ference Series, Vol. 351, Astronomical Data Analysis Soft¬ 
ware and Systems XV, Gabriel C., Arviset C., Ponz D., 
Enrique S., eds., p. 112 
Berlin E., Arnouts S., 1996, A&AS, 117, 393 
Berlin E., Mellier Y., Radovich M., Missonnier G., Didelon 
P., Morin B., 2002, in Astronomical Society of the Pacific 
Conference Series, Vol. 281, Astronomical Data Analysis 
Software and Systems XI, Bohlender D. A., Durand D., 
Handley T. H., eds., p. 228 

000 , 


Botticella M. T., Smart! S. J., Kennicutt R. C., Cappellaro 
E., Sereno M., Lee J. C., 2012, A&A, 537, A132 
Castelli F., Kurucz R. L., 2004, astro-ph/0405087 
Chini R., Hoffmeister V. H., Nasseri A., Stahl O., Zinnecker 
H., 2012, MNRAS, 424, 1925 
Dalcanton J. J. et ah, 2012, ApJ, 200, 18 
Davidson K., 1987, ApJ, 317, 760 
Davidson K., Humphreys R. M., 1997, ARA&A, 35, 1 
Dessart L., Livne E., Waldman R., 2010, MNRAS, 405, 
2113 

Dolphin A. E., 2000, PASP, 112, 1383 
Draine B. T., 2011, Physics of the Interstellar and Inter- 
galactic Medium 

Draine B. T., Lee H. M., 1984, ApJ, 285, 89 
Elitzur M., Ivezic Z., 2001, MNRAS, 327, 403 
Foellmi C. et ah, 2008, Rev. Mexicana Astron. Astrofis., 
44, 3 

Freedman W. L. et al., 2001, ApJ, 553, 47 
Gerke J. R., Kochanek C. S., Stanek K. Z., 2015, MNRAS, 
450, 3289 

Horiuchi S., Beacom J. F., Kochanek C. S., Prieto J. L., 
Stanek K. Z., Thompson T. A., 2011, ApJ, 738, 154 

000-000 


0000 

















Imposter SN 1997hs 15 


Humphreys R. M., Davidson K., 1994, PASP, 106, 1025 
Humphreys R. M., Davidson K., Smith N., 1999, PASP, 
111, 1124 

Ivezic Z., Elitzur M., 1997, MNRAS, 287, 799 
Ivezic Z., Nenkova M., Elitznr M., 1999, astro-ph/9910475 
Khan R., Kochanek C. S., Stanek K. Z., Gerke J., 2015, 
ApJ, 799, 187 

Khan R., Stanek K. Z., Kochanek C. S., 2013, ApJ, 767, 
52 

Kochanek C. S., 2009, ApJ, 707, 1578 
Kochanek C. S., 2011a, ApJ, 741, 37 
Kochanek C. S., 2011b, ApJ, 743, 73 
Kochanek C. S., 2014a, MNRAS, 444, 2043 
Kochanek C. S., 2014b, arXiv: 1407.7856 
Kochanek C. S., 2014c, ApJ, 785, 28 
Kochanek C. S., 2015, MNRAS, 446, 1213 
Kochanek C. S., Beacom J. F., Kistler M. D., Prieto J. L., 
Stanek K. Z., Thompson T. A., Yiiksel H., 2008, ApJ, 
684, 1336 

Kochanek C. S., Szczygiel D. M., Stanek K. Z., 2011, ApJ, 
737, 76 

Kochanek C. S., Szczygiei D. M., Stanek K. Z., 2012, ApJ, 
758, 142 

Krist J., 1995, in Astronomical Society of the Pacific Con¬ 
ference Series, Vol. 77, Astronomical Data Analysis Soft¬ 
ware and Systems IV, Shaw R. A., Payne H. E., Hayes 
J. J. E., eds., p. 349 

Krist J. E., Hook R. N., Stoehr F., 2011, in Society of 
Photo-Optical Instrumentation Engineers (SPIE) Confer¬ 
ence Series, Vol. 8127, Society of Photo-Optical Instrn- 
mentation Engineers (SPIE) Conference Series 
Knshnir D., Katz B., 2014, arXiv: 1412.1096 
Earners H. J. G. L. M., de Groot M., Cassatella A., 1983, 
A&A, 123, L8 

Li W., Filippenko A. V., Van Dyk S. D., Hu J., Qiu Y., 
Modjaz M., Leonard D. C., 2002, PASP, 114, 403 
Mathis J. S., Rumpl W., Nordsieck K. H., 1977, ApJ, 217, 
425 

Pan K.-C., Ricker P. M., Taam R. E., 2013, ApJ, 773, 49 
Pastorello A. et ah, 2004, MNRAS, 347, 74 
Pejcha O., Thompson T. A., 2015, ApJ, 801, 90 
Prieto J. L. et ah, 2008, ApJ, 681, L9 
Puls J., Vink J. S., Najarro F., 2008, A&A Rev., 16, 209 
Saha A., Sandage A., Tammann G. A., Labhardt L., Mac- 
chetto F. D., Panagia N., 1999, ApJ, 522, 802 
Schlafly E. F., Finkbeiner D. P., 2011, ApJ, 737, 103 
Schlegel D. J., Finkbeiner D. P., Davis M., 1998, ApJ, 500, 
525 

Shappee B. J., Kochanek C. S., Stanek K. Z., 2013, ApJ, 
765, 150 

Smart! S. J., 2009, ARA&A, 47, 63 
Smith N., Humphreys R. M., Gehrz R. D., 2001, PASP, 
113, 692 

Smith N., Li W., Silverman J. M., Ganeshalingam M., Fil¬ 
ippenko A. V., 2011, MNRAS, 415, 773 
Smith N., Owocki S. P., 2006, ApJ, 645, L45 
Spiro S. et ah, 2014, MNRAS, 439, 2873 
Thompson T. A., Prieto J. L., Stanek K. Z., Kistler M. D., 
Beacom J. F., Kochanek C. S., 2009, ApJ, 705, 1364 
Treffers R. R., Peng C. Y., Filippenko A. V., Richmond 
M. W., Barth A. J., Gilbert A. M., 1997, lAU Circ., 6627, 
1 


Turatto M. et ah, 1998, ApJ, 498, L129 
Ugliano M., Janka H.-T., Marek A., Arcones A., 2012, ApJ, 
757, 69 

Van Dyk S. D., Filippenko A. V., Chornock R., Li W., 
Challis P. M., 2005, PASP, 117, 553 
Van Dyk S. D., Matheson T., 2012, in Astrophysics and 
Space Science Library, Vol. 384, Astrophysics and Space 
Science Library, Davidson K., Humphreys R. M., eds., p. 
249 

Van Dyk S. D., Peng C. Y., Barth A. J., Filippenko A. V., 
1999, AJ, 118, 2331 

Van Dyk S. D., Peng C. Y., King J. Y., Filippenko A. V., 
Treffers R. R., Li W., Richmond M. W., 2000, PASP, 112, 
1532 

Wagner R. M. et ah, 2004, PASP, 116, 326 
Weis K., Bomans D. J., 2005, A&A, 429, L13 
Woosley S. E., Heger A., 2015, in Astrophysics and Space 
Science Library, Vol. 412, Astrophysics and Space Science 
Library, Vink J. S., ed., p. 199 
Woosley S. E., Weaver T. A., 1995, ApJS, 101, 181 


000 , 000-000 



